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1 Introduction 
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j^ , Let An — C(a;i, . . . , a;„, 9i, . . . , dn) be the Weyl algebra of order n over the com- 

plex numbers C and C[d] — C[di, . . . , dn] the subring of An of linear differential 
operators with constant coefficients. 

Let A = {aij) be an integer d x n- matrix of rank d. We denote by I a C C[d] 

^^^ , the toric ideal associated to A: i.e. I a is the ideal generated by the set 

f~^ ■ where Q'^ means "transpose". 

r^ [ We denote by 9 the vector {9i, . . . , 9n)^ with 9i — Xidi. For a given /3 = 

(/3i, . . . ,/3d)"^ £ C'' we consider the colmiin vector (in Af^) A9 — (i and we denote 
^^ \ by {A9 — /?} the left ideal of An generated by the entries of A9 — j3. 

f^ ■ Following Gel'fand, Kapranov and Zelevinsky Q, we denote by Ha{0) the 

"jrt ! left ideal of An generated by I a U {A9 — /3) . It is called the GKZ-hypergeometric 

system associated to the pair (A, /3). The quotient TiA{P) = An/HA{[i) is a 
holonomic A„-module (see e.g. |17|). 

If the toric ideal I a is homogeneous, i.e., if the Q-row span of A contains 

k><( \ (1, . . . , 1), it is known (Q, see also [|l^) that HAiP) is regular holonomic and 

j_j ■ the book |l^ is devoted to an algorithmic study of such systems. Especially, the 

Cu I book gives an algorithmic method to construct series solutions around singular 

points of the system. 

In this article, we would start a study of singularities of GKZ-hypergeometric 
systems for non homogeneous toric ideals Ia by treating the "first" case when 
d = 1, A = (ai,a2, . - . ,a„) £ Z", ai = 1. We evaluate the geometric slopes 
of TiAiP) by successive restrictions of the number of variables. The slopes 
characterize Gevrey class solutions around singular locus. 

Our evaluation of the slopes is done as follows: (1) We translate Lau- 



rent and Mebkhout's theorem |12 on restrictions and slopes of P-modules 
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into an algorithm to evaluate the slopes by utilizing the results [^ and | p^ . 
(2) Apply our general algorithm to the hypergeometric system associated to 
^ = (1, 02, . . . , a„). This system has a lot of nice properties and our algorithm 
outputs the slopes without computation on computers. 

2 Micro-characteristic varieties 

In this section following Laurent ||l^ , we describe micro-characteristic varieties 
for a given X'-module. We will state a result of Laurent and Mebkhout |^, 
CoroUaire 2.2.9 ] (see also (@, p. 125) and (jlll, p. 42)), allowing to reduce 
our general problem of evaluating the slopes to less variables. 

In this section X = C" and Vx = I? is the sheaf of linear differential 
operators with holomorphic function coefficients. Let M be a coherent "Dx- 
module. Recall that the characteristic variety of M (denoted by Ch(Af )) is a 
analytic sub-variety of the cotangent bundle T*X. 

Suppose that y C X is a smooth hypcrsurface. We say that Y is non- 
characteristic for M if TpX n Ch(M) C TJX. Here T^X is the conormal 
bundle to F in X and T^X is the zero section of T*X. 

Now, following Laurent B, Efl, we shall define the notion of non-micro- 
characteristic variety for M . To simplify the presentation we will assume that 
(xi, a;2, . . . , x„) are local coordinates in X and Y is defined by x„ — 0. We 
denote by (xi, . . . , a;„, ^i, . . . , ^„) local coordinates in T*X. Sometimes it will 
be useful to write xi = yi, . . . , a;„_i = 2/„_i, a::„ = i, ^i = 771, . . . , ^„_i = ?7„_i 
and S,n = T. 

Let us denote by A the conormal bundle of Y in T*X (i.e. A = TpX). So, in 
local coordinates A = {{y,t,ri,T) G T*X \t = t] = 0}. Here, rj = {iji, . . . ,rjn-i) 
and y — (yi, . . . , j/„_i). We denote by {y,T,y*,T*) local coordinates on the 
cotangent bundle T*A. 

We denote by V,(2?) (or simply by V) the Malgrange-Kashiwara filtration 
associated to F on D and by Ft{'D) (or simply by F) the order filtration on 
v. For a given rational number p/q > we denote by Lp/g the filtration on 
V defined by pF -\- qV. The Lp/^-order of a monomial y^t^df^d^ is equal to 
p{\f3\-\-k)-\-q{k — l), where |/3| = /3i 4- • • ■+/3„_i. We will simply write L = L^jq 
if no confusion arises. 

For p > the associated graded ring gr^(X') is canonically isomorphic to 
T^*0[T*A\ [|lO[ P- 407], where n : T* K — ^^ A is the canonical projection and 
0[T'A] denotes holomorphic functions on T* K which are polynomials on the 
fibers of TT. In local coordinates, gr^(X'o) is expressed as C{y}[T, y*,T*] where 
Vq is the stalk of T> at the origin. 

Given a differential operator 

P^Y. Po^iPky^t'dld'l 

aipk 

the L-order of P is the maximun value oip{\j3\+k) + q{k~l) over the monomials 



of p. For p > we define the L-principal symbol of P by 

where the sum is taken over monomials with maximal L-order. The i-principal 
symbol of P is an element of gr^(P) and then is a fmiction on r*A. In the 
classical case, i.e., for L — F, gr^(I?) is identified with C{a::}[^i, . . . ,^„] = 
C{y,t}[r],T] and the i^-principal symbol of P is simply denoted by a^{P) = 
^Pai^kV^t^'n^ T^ where the sum is taken for |/3| + k maximum. 

For each left ideal I C V we denote by a^ (/) the ideal of gr-^ {V) generated 
by the set of cr^(P) for Pel. 

For each L-filtration on T> we associate a "good" L-filtration on M, by means 
of a finite presentation. The associated gr-^ (I?)-modulc gr^ (M) is coherent (see 
[ pO| 3.2.2]). The radical of the annihilating ideal Anngrt(x)jf)(gr'^(Af)), which 
is independent of the "good" filtration on M, defines an analytic sub- variety of 
T*A. This variety is called the L- characteristic variety of M and it is denoted 
by Ch^(M). 

Suppose now Z C X is a smooth hypersurface transverse to Y. Suppose 
for simplicity Z is defined in local coordinates by j/i =0. The conormal space 
A' := Ty^z^ i^ ^ smooth subvariety of A = TyX defined in local coordinates 
by yi = 0. So T^,A is the subvariety of r*A defined in local coordinates by 
2/1 = y'* = r* = 0, where y' = {y2, ■ ■ . ,2/„-i). 

Definition 2.1 [||, iQ We say that Z is non-micro-characteristic of type L for 
M if T^/A n Ch^(M) is contained in T^A. Sometimes we will say that, if this 
condition holds, Z is non L -micro- characteristic for M. 

The sheaf of rings gr^(2?) is endowed with two graduations: First one is 
induced by the F-filtration and the second one is induced by F-filtration. Recall 
Laurent's definition of slope of a coherent P-module M . 

Definition 2.2 [nOJ The rational number —p/q is said to be a slope of M w.r.t. 
Y at the origin if and only if the radical of the ideal Anngi.t(i5^)(gr^(M)) is not 
bi-honiogeneous for F nor V graduations. 

An important consequence of the work ||ll| (see Th. 2.4.2) is what follows: An 
holonomic Px-module M is regular w.r.t. Y at the origin if and only if M has 
no slope w.r.t. Y at the origin. We will use this fact freely in the text. 
Remark: In [Q, oo and {F and V) are included in the set of the slopes. We 
do not include them in the set of the slopes in this paper. 

Finally the following result by Laurent and Mebkhout allows induction on 
the number of variables to calculate slopes. 

Theorem 2.3 |12, CoroUaire 2.2.9] Let M he a holonomic Vx -module. Let Z 
and Y be transverse smooth hypersurfaces on X such that Z is non L-micro- 
characteristic (for all p > 0,q > 0) for M . Then the slopes of M w.r.t. Y equals 
the slopes of M' w.r.t. Y' ^ Z ClY where M' is the restriction of M to Z. 



This is a deep result in I'-inodule theory. Its proof uses the algebraic-analytic 
comparison theorem (pil Th. 2.4.2 ]) and a Cauchy-Kowalewska theorem for 
Gevrey functions w.r.t. Z (|l|, CoroUaire 2.2.4]; see also Q Th. 6.3.4].) 

3 Computing slopes by reducing the number of 
variables 

We have introduced the notion of the slopes and the invariance of them under 
restrictions satisfying a condition on L-characteristic varieties. In the sequel, 
we assume that our ideal is that of the Weyl algebra A^. Constructions in 
sheaves such as restrictions, micro-characteristic varieties in the previous section 
can be done via constructions in the Weyl algebra as we usually see in the 
computational I?-module theory. 

We are interested in computation of the slopes. The slopes of A^jl (at the 
origin) along a:„ = can be computed by the ACG algorithm introduced in 
[0. In this section, translating Laurent and Mebkhout's result into computer 
algebra algorithms, we will give a preprocessing method for the ACG algorithm 
to accelerate the original. The preprocessing is useful for a class of inputs 
including GKZ hypergeometric ideals as we will see in Section H. Let us firstly 
recall the ACG algorithm. 

A weight vector is an element W = (ui, . . . , u„, wi, . . . , f„) G R^" such that 
Ui-\- Vi>Q for all i. This weight vector W — {u,v) induces a natural filtration 
on An and it is called the W^-filtration. The associated graded ring is denoted by 
gT^{An) and for each left ideal I C An the associated graded ideal is denoted 
by inw{I) or in(„ „•)(/). Here, the initial ideal iii(^u.v){I) is the ideal generated by 
in(.^ „)(/), / G J in gr^(y4„). When Ui + Vi > 0, it is an ideal in the polynomial 
ring of 2n variables: gr^(A„) = C[a:i, . . . ,Xn,£,i, ■ ■ ■ ,^n]- The initial ideal of 
/ with respect to the weight {u, v) is generated by the {u, i;)-initial terms of a 
Grobner basis of / by an order which refines the partial order defined by {u, v). 
See, e.g., 0, Theorem 1.1.6]. 

Consider the filtration L — pF + qV,p > 0,q > introduced in the previ- 
ous section. The ideal cr^{I), which gives the L-characteristic variety, can be 
expressed in terms of the initial ideal as follows: 

a^{I) = gr^(P) • iuf (/)UiH^yi,...,x,._iH^y„_i,x„H^r*,6H-!yi*,...,?„-iH-y:_i,C„H^-r, 



£ = p(0,...,0, !,...,!) +(7(0,..., 0,-1,0,..., 0,1). 

For two weight vectors W and W' and a term order <, we denote by <w,W' 
the order 

x"d^ <w,w' x^d^ 
^ W ■{a,(3)<W -{a.h) 

or W ■{a,l3)^W ■ (a, h) and W ■ (a, P) <W' ■ (a, b) 



■ W* ■ (a, I3)^W* ■ (a, b) for both W* = W, W and x'^d'^ < x'^d'' 



To each differential operator P — ^PafiX°'d^ G An we associate tlie Newton 
polygon N{P) of P (w.r.t. Xn — 0) defined as the convex hull of the subset of 

U (|/3|,/3„-a„) + (-N)2. 



Let / be a left ideal in A„. As we said in Definition 2.S the notion of slope 
of a differential system was introduced by Y. Laurent |lC|. Let us give here 
a slight different but equivalent definition: the number r, — cxd < r < Oisa 
geometric slope of / (or of A„/J) w.r.t. a;„ = if and only if ^cr(^'')^+^(/) 
is not bihomogeneous with respect to the weight vectors F = (0, . . . , 0, 1, . . . , 1) 
and V = (0, . . . ,0, — 1,0, . . .,0, 1). Following 0, we say that the number r, 
— oo < r < is an algebraic slope of I (or of An/ 1) if and only if cr(~'')^+^(/) is 
not bihomogeneous with respect to the weight vectors F and V . The geometric 
slope is simply called the slope in this paper if confusion does not arise. Note 
that we may consider in^ (/) instead of a^ (/) as far as we are concerning about 
homogeneity. For algebraic or geometric slope r, the weight vector L — (— r)F + 
V lies on a face of the Grobner fan of / ( 0] , [O ) , which yields the following 
algorithm. 

Algorithm 3.1 (|], ACG algorithm) 

Input: G — {Pi, . . . , Pni} (generators of an ideal /) 

Output: All algebraic and geometric slopes of An/ 1 w.r.t. x„ = at the origin. 

geometric_slope — 0; algebraic_slope — 0; 

F=(0,.. .,0,1,. ..,1);V=(0,.. .,0,-1,0,. ..,0,1); 

p = 1; q = 0; slope = —cxd; previous_slope = slope; 

while (slope! = 0) { 

L = pF + qV; 

G = a Grobner basis of / with respect to the order <l,v\ 

slope = the minimum of and 

{ the slopes r of the Newton polygon N{P) | P £ G, r > previous_slope} 

if slope = , then return(algebraic_slope and geometric_slope). 

if cr^{G) is not homogeneous for F nor V then { 
algebraic_slope = algebraic_slope U {slope} 

} 

if yoHfi) is not homogeneous for F nor V then { 
geometric_slope = geometric_slope U {slope} 

} 

p = numerator(|slope|); q = denominator(|slope|); 

previous_slope = slope; 

} 

return (algebraic_slope and geometric_slope) 



Here, we use the convention F :— ooF+V . Note that the ideal J ^ (/i , . . . , /,„) C 
C[a:i, . . . , a;„, ^1, . . . , ^„] is homogeneous for the weight [u, v) G Z^" if and only 
if all (m, u)-homogeneous subsums of /,; belong to the ideal J. 



For ordinary differential equations, tfie ACG algoritlim is notliing but tlie well- 
known Newton polygon method; see two examples below. 

Example 3.2 Put G = {xf+^^i +p}, n = 1, p G N. Then, the ACG algorithm 
returns {—p}- ( exp(a;7'') is a classical solution of G.) 

Example 3.3 Put G = {2a;i(2;iai)2 + xldi + 1}, n ^ 1. Then, the ACG 

— 1/2 

algorithm returns { — 1/2}. ( exp{x^ ) is a classical solution of G.) 

Next example of two variables is generated by a computer algebra system kan/kO 



Example 3.4 We consider the GKZ hypergeometric ideal / associated to the 
matrix A = (1,3) and f3 — —3. We will compute the slopes of A2/I at the origin 
along a;2 = by the ACG algorithm. The Grobner basis of / for the weight 



Xi X2 di 



1 



IS 



[ xl*Dxl+3*x2*Dx2+3 , -Dxl~3+Dx2 , -3*x2*Dxl"2*Dx2-5*Dxl"2-xl*Dx2 , 
-9*x2"2*Dxl*Dx2"2-36*x2*Dxl*Dx2+xl"2*Dx2-20*Dxl , 
-27*x2-3*Dx2-3-189*x2-2*Dx2-2-xl-3*Dx2-276*x2*Dx2-60 ] 

Here, Dxi and xi stand for di and Xi respectively. The Newton polygons 
N{Pys are figured below. 



^ T 



=P^ 



=Pt 



=F^ 



— II m • ^ 



From the Newton polygons, —1/2 is the candidate of the first slope. Next, we 
compute the Grobner basis for the weight (0, —2, 1, 3) = (0, 0, 1, l)-|-2(0, —1, 0, 1). 
The Grobner basis is 



[ xl*Dxl+3*x2*Dx2+3 , Dx2-Dxl"3 ] 

The radical is generated by 

[ -xl*Dxl-3*x2*Dx2 , -Dxl"3+Dx2 , 3*x2*Dxl"2*Dx2+xl*Dx2 , 
-9*x2"2*Dxl*Dx2"2+xl"2*Dx2 , 27*x2"3*Dx2"3+xl"3*Dx2 ] 

Here, xl = yi, x2 = r*, Dxl = yl, Dx2 == — r. It is not bi-homogeneous and 
then —1/2 is a geometric and algebraic slope. 

By looking at the two Newton polygons of xl*Dxl+3*x2*Dx2+3, Dx2-Dxl~3, 
we see that there is no more slope that is larger than —1/2. Then, the ACG 
algorithm terminates here. 

The ACG algorithm requires a repetition of Grobner basis computations in 
the Weyl algebra of 2n variables to evaluate the slopes of An/ 1. However, if 
Xi = is non-micro-characteristic of type L = {—r)F + V for all — oo < r < 
and the restriction of An /I to Xi ^ is singly generated, we can preprocess the 
input so that the input ideal for the ACG algorithm lies in the Weyl algebra of 
2(n — 1) variables. The correctness of the following algorithm can be shown by 



Laurent and Mebkhout's theorem 2.3 



Algorithm 3.5 (Computing slopes with a preprocessing) 

Step 1: check if Xj = is non- micro-characteristic of An /I for all types L by 



calling Algorithm 3.6 



Step 2: Compute the restriction of An /I to Xi = and check if it is expressed 

as D' / 1' where D' is the Weyl algebra of 2{n — 1) variables. 

Step 3: If we failed either in Step 1 or in Step 2, then apply the ACG algorithm 

for /. 

If we succeeded both in step 1 and in step 2, then try to reduce more variables 

or apply the ACG algortihm for /'. 

We can compute restrictions of a given 2?-module by using Oaku's algorithm 
[y_5|. This algorithm is implemented in computer algebraic systems Macaulay2 
and Kan |Q, iQ, |Q. Therefore, the remaining algorithmic question for the 
preprocessing is to determine the range of type L = {—r)F -\- V for which 
Xi = {i < n — 1) is non-micro-charactcristic. It follows from the definition of 
non-micro-characteristic that the question is nothing but to find the segment 
(— oo,ri) such that 

V(a(-'-)^+^(/), yl . . . ,y*_i,2/„ y*+i, . . . ,y,:_i,T*) C Viyl . . . , y*n.„T*) = T^A 

for r € (— oo, ri). Here, V(/i, . . . , fm) is the affine variety defined by the poly- 
nomials /i, . . . , fni- The inclusion condition can algebraically rephrased as 



yin(_r)F+v(-^),Cii ■ ■ ■ :Cj-i,a;i,^i+i, ■ ■ ■ ,Cn ^ ^i- 

Since the Grobner fan is a finite union of Grobner cones H, the range can be 
determined by a similar method with the ACG algorithm. 



Algorithm 3.6 

range_of _nonMC(_ff, tq) 

Input: H is a finite set in An, tq is a negative number or — oo 

Output: ri such that xi = is non-niicro-characteristic of type {—r)F + V 

for An/ An ■ {H], for r e [ro,ri) . 

previous_slope = slope = tq; 

F=(o,. .. ,0,1, •••,i);v=(o,... ,0,-1,0,. ..,0,1); 

while (slope! = 0) { 

p = nunierator(|slope|); q ~ denominator(|slope|); 
L = pF + qV; 

G = a Grobner basis of H with respect to <l\ 
if ^/ {<^''{H),yi,yl, . . . ,yl_^,T*) B y* and 



(^pF+(q+s)v(^)^ yi,2/2% ■ ■ -^Vn-i^r*) 3 yl 
then { 

previous_slope = slope; 
slope — the minimum of and 

{ the slopes r of the Newton polygon N{P) | P G G, r > previous_slope} 
} else { 

rcturn(slope); 

} 
} 

return(O); 

Here, e is a sufficiently small positive rational number so that pF + (9 + £)V 
lies in the interior of a Grobner cone. 

When ro = —00, we use the convention numerator(|ro|) = 1 anddenominator(|ro|) = 
and P-non-micro-characteristic means that it is non-characteristic in the clas- 
sical sense. 

Example 3.7 Suppose n > 2. range_of .nonMCd^j' — 82}, —00) returns 0. 

If the function range_of _nonMC(/, —00) returns 0, then xi = is non-micro- 
characteristic of type pL + qV {p > 0, q > 0) for An/I- We note that it is not 
always a clever strategy to call the function with the full set of generators. In 
fact, ii xi = is non-micro-characteristic of type L for An/ J, then it is non- 
micro-characteristic of type L for An /I for any I D J. Therefore, for step 1, it 
is sometimes more efhcient to call the function range.of _nonMC for a subset of 
the generators of the input ideal as we will see in case of GKZ hypergeometric 
ideals in the next section. 

4 Computing slopes of ^^(las, ,«„)(/?) 

Put A = (1, 02, • • • , a„), ai = 1 < 02 < • ■ • < a„. We will evaluate the slopes 
of the GKZ hypergeometric ideal Ha (/3) associated to the 1 x n matrix A and 
/? G C by using the general algorithm given in Section 0. To apply this al- 
gorithm, we need to find non-microcharacteristic varieties and compute the 



restrictions of TtA{P) to these varieties. When /i, . . . , /„ are polynomials in 
C[a:i, . . . ,a:„,^i, . . . ,^„], we denote by V(/i,...,/m) the affine subvaricty in 
C2" defined by the f,. 

The following theorem can be shown by a standard method of Koszul com- 
plex (0, §). 

Theorem 4.1 The characteristic variety ofTiAiP) ^^ ^(fii ■ • ■ jCn-ij^^nCn)- ^'^ 
particular, the singular locus ofTiAiP) ^^ ^n — 0. 

Note that there is no slope along Xi = 0, l<i<r7, — 1, which can be shown 
easily. 

Recall that i^ = (0, . . . ,0, 1, • • ■ , 1) and X^ = (0, . . . ,0, -1,0, . . . ,0, 1). For 
positive number p and a non-negative number q, we define a weight vector 
L = pF + qV. In section g, we explained the notion of non-micro-characteristic. 
When the variety is yi = Xi — 0, this notion is rephrased as follows: for a 
given left A„-module An/I, the hyperplane j/^ = (1 < z < ti — 1) is called 
non-micro-characteristic of type L when 



^{<T'^{I),y^,y*,{J^^),r*)^y* 



Proposition 4.2 For the hypergeometric A„-module Ti.A{P), the variety yi = 0, 
(1 < i < n — 2) is non-micro- characteristic of type L for all L = pF-\-qV, p > 0. 

Proof: Consider 5f^ — d"^' G Ha{P). For all L and for i < j < n — 1, we have 
a^{d^' — d":^) — {y*Y^ , which implies that y^ = is non- micro-characteristic of 
type L. 

Now, let us apply the second step of the algorithm to evaluate the slopes, i.e., 
we will compute the restriction of 7iyi(/3) to y^ = x^ = 0, (1 < z < n — 2). Let 
s be an indeterminate. Consider the ideal Ha[s\ in A„[s] generated by A6 — s 
and Ia- 

Theorem 4.3 We have a left D'[s]-module isomorphism 

Ar,[s]/iAn[s]HA[s]+x,An[s])-D'[s]/D'[s]HA'[s], i^l (1) 

Here, D' = C(xi, . . . , x,_i, x^+i, . . . , a;„, i9i, . . . , (9i_i, 9^+1, . . . , 9„) and A' = 
(1, a2, . . . , ai_i, fli+i, . . . , a„). 

Proof: Fix an order -< such that di >- dn >- 9„-i )^ ■ ■ ■ y di. Then, 

dn-d'^-, dn-i-d'i-\ •••,92-a?^ 

is the reduced Grobner basis of I a with respect to y. 

For each i {2 < i < n), hy applying the same method with the proof of 
[ p^Tt Th 3.1.3], we can prove that in(_g. e.')(iJ^(s)) is generated by A6 — s and 
ine,(9j -9i'), (j = 2, ...,n). 



Define the 1/^-filtration Fk[s] of An[s] by 

PkM = {Y. a^p^x^d^s^ I (/3 - a) • e, < fc} . 

We compute the restriction of A„[s]/i/^(s) to a;,; = by Oaku's algorithm (see, 
e.g., ||l^. Theorem 5.2.6, Alg. 5.2.8]). Since inj.g. e.-)(a;i9i — Xidi') — Xidi, the 
6-function oi Ha{s) along Xi = is b{p) — p. Therefore, by |]1^, Theorem 5.2.6], 
we have the following isomorphism as left I?'[s]-modules: 

An[s]/{HA{s)+X,An[s]) 

n 

~ Fo[s]/{Fo[s]{Ae -s)+ J2 ^oW(<9, - d^,^) + F_,[s]{d, - d^^) + x,F^[s]) 
^ D'[s]/lD'[s]{Ae-s)+ J2 D'i^M-d'in 

D 

We can specialize s to any complex number /3. So, we have 

A„/iArMAW + X,An) ^ D'/D'Ha'{P). 

which means that the restriction of TCa (/3) to x^ = can be exactly expressed 
in terms of the GKZ system for smaller A. By applying our algorithm ^.5| of 
computing slopes by reduction of number of variables to variables X2, ■ ■ ■ , Xn-2, 



we obtain the following theorem by Proposition 4.2 and Theorem W^ 



Theorem 4.4 The geometric slopes oJ1-La{P) along x^ = at the origin and 
7i(x.a„_i.a„)(/3) along x^ = at the origin coincide. 

Example 4.5 We note that (the algebraic slopes) ^ (the geometric slopes) in 
general. For example, let us apply the ACG algorithm to get the algebraic 
slopes of i?A(— 30) for A = (1, 3, 7). This ideal is generated by 

a;i9i + 3x292 + Tx^da + 30, df - 82, ^dfd2 + 83, 9| - dfda. 

The output is 

{-1,-3/4,-1/2}. 

On the other hand, if we apply the ACG algorithm to get the geometric slopes, 
the output is {—3/4}. 

Example 4.6 (the slopes oi H^i,a„_i,a„){f^)) ¥" (the slopes of 'W(i^a^)(/3)) in gen- 
eral. 

Let us take the example: A — (1,3,7). Consider the hypergeometric As-module 
A3/I where / — i?(i 3 7)(— 30). As we have seen, the slope of this system along 
X3 = is {-3/4}. 
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Consider iriL (/) ior L = F + 4V. By computing Grobncr basis with respect 
to L, we can see that 

V(ini(/))=V(a;2,ei,e3)UV(Ci,6,6)- 

It is not included in V(^2) Hence X2 = is micro-characteristic of type L and 
we cannot apply for the restriction criterion. 

The condition "non-microcharacteristic for all the filtration pi^ + gV^" cannot 
be taken a way to evaluate the slopes by the restriction. In fact, it can be easily 
checked by the ACG algorithm that the set of the geometric slopes oi'H(^i^a„)iP) 
is equal to {1/(1 — a„)}. Hence, the set of the slopes of 7i(i,7)(/3) is { — 1/6}, 
which is not equal to {—3/4}. 

We have shown that the computation of the slopes of Ha (/3) is reduced to 
the three variables case. The slopes of this case are as follows. 

Theorem 4.7 

(the slopes of H(i,a„_i,a„)(/3)) = {a„_i/(a„_i - a„)} 
Proof: We fix some notation: 

1. A = (1, a, b) e Z^ and 1 < a < 6. 

2. Pi = dl -d2,P2= d'l - ds, Pg = dl - d^, Pi = xidi + ax2d2 + bx^da - /3 

3. Let A be the linear form with slope —a/{b — a) (i.e. L ~ aF + {b — a)V). 

4. Let L, L' be linear forms. We say that L > L' \i slope(P) > slope(L'). 

5. We will write yi — xi, j/2 — a^2, t = a^3- 

The operators Pi, Pj, P3, P4 are in iJ = Ha{(3)- Let A be a„_i/(a„ — a„_i). 
Then, we have the following claims. 

1. For all linear form L we have ct^(Pi) ~ (?7i)" and so 77* G \Ja^{H) for all 
L. 

2. For all linear form L we have a^{P4) — yirj\ + ay2'72 + bT*{~T) 

3. For all linear form L > A we have cr^(P3) = (r)" and so r £ \/(7^{H) for 
all L > A. 



4. So, for all i > A we have Ch^(H) C T*^^qC^\JT^sC^ and then ^JaHH) 
is bi- homogeneous and L is not a geometric slope of Ti. 

5. On the other hand we have, for P < A, cr^(P3) = {flt)^ ■ Then 773 G 
^(7^{H) and Ch^(H) C T^^^C^ U T^:,C^ . So L is not a geometric slope 
of H. 

6. So, the only possible geometric slope of Ti is A. 
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Now, suppose that A is not a slope. Then, there is no slope, which implies 
that the L-characteristic variety Ch {T-Ca{0)) is the same for all L = pF + qV, 
p, g > by H and 0, Th 3.4.1]. It follows from 

[PuP2] = 0, [Pi, Pi] = aPi, [P2, Pi] = bP2 

and the Buchberger algorithm that {Pi, P2, P4} is a Grobner basis for the order 
defined by the weight vector L = (0, 0, —N, 1, 1, N+1), N > b and a tie-breaking 
term order such that X2 ^ 83 y di y 82- Therefore, the initial ideal ini(iJ^(/?)) 
is generated by ^°, ^3, a;i^i-|-aa;2^2 + ^2;3^3 and hence the L-characteristic variety 
is equal to T*^qC^ U T^^C^. This fact contradicts that the L-characteristic 
variety is the same for all L. |] 

Finally, let us remark on an analytic meaning of slopes. Let X — C" and 
Y — {x\xn = 0} C X. We denote by Orrpr; = C{a;i, . . . ,x„_i}[[a:„]] the formal 

completion of Ox along Y. For each real number s e [l,-|-oo) we denote by 
Ox\y{s) the subsheaf of 0-— - of Gevrey functions of order s (along Y). The 

sheaf 0x1^(1) is the restriction Ox\y and, by definition, we write Ox\y{+oo) = 
O-— -. For any holonomic 2?x-niodule M, Mebkhout associates in ||lj] the 

sheaf Irry(s) (M) as the solution shea,i'RHomj){M,Ox\Y{s)/Ox\Y)- One fun- 
damental result in the irregularity of 2?- modules is the fact that Irry (s)(M) is 
a perverse sheaf, for any s (see |l3]). These sheaves define a filtration of the 
irregularity of M along Y, i.e. Irry(M) := Irry(-|-oo)(M). The main result of 
[ pr] is that 1/(1 — s) is a slope of M w.r.t. Y if and only if s is a gap of the 
graduation defined by the filtration on the irregularity. In other words, 1/(1 — s) 
is a slope if and only if Irry(s)(Af)/IrrY(< s){M) ^ 0. 

5 Rational Solutions and Reducibility 

Our ultimate aim of studying the slopes of HA{f3) is to get a better under- 
standing on solutions of this system. We are far from the goal, but to this end, 
it will be useful to present some facts on classical solutions and a relation to 
generalized confluent hypergeometric functions. 

In case of the hypergeometric ideal associated to homogeneous monomial 
curves, Cattani, D'Andrea, Dickenstein M studied rational solutions and re- 
ducibility of the system. We will study rational solutions and reducibility of our 
system. 

Theorem 5.1 Any rational solution of the hypergeometric system -ff(i,a2,...,a„)(/^) 
is a polynomial. It has a polynomial solution if and only if f3 € N — {0, 1, . . . , }. 
The polynomial solution is the residue o/exp (^ Xit""^) t^^ at the origin t ~ 0; 

X»p(i:-'-)'-'t 

Here, C is a circle that encircles the origin in the positive direction. 
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Proof: Since the singular locus of Ha{P), {A = (ai, 02, . . . , a„), ai — 1) is 
Xn — 0, any rational solution / is a Laurent polynomial with poles on a;„ = 0. 
Take a weight vector w = (0, 1, 1, . . . , 1). Then, we have inu,(JA) = (92, . . . , i9„). 
The initial term iu-wif) is annihilated by 

'^aiO^-P, d2,...,dn. 

Therefore, x^ — iuwif)- This implies /3 G N, because / has a pole only on 
Xn = 0. 

Take a Z-basis of Ker (Z" 4 Z) as 

(-a2, 1,0,..., 0), (-as, 0, 1, 0, . . . , 0), . . . , (-a„, 0, . . . , 0, 1) 

to construct series solutions. Since [17| Prop. 3.4.1] holds for non-homogeneous 
A as well, formal series solution g of Ha{P) satisfying in^(5) = x^ can be 
uniquely expressed as 

When /3 e N, it is a polynomial. The rest part of the theorem is easy to show. 



Let R be the ring of differential operators of n variables with rational function 
coefhcients over k = C. A left ideal J of i? is called irreducible when J is a 
maximal ideal in R. We will study the reducibility of R ■ Ha (/3) . 

We assume that J is zero-dimensional, i.e., r = dimc(i:) R/ J < +cxd. Let 
V = V{J) be the vector space of holomorphic solutions of J on a simply con- 
nected open set contained in the non-sigular domain of J. It is known that 
dime V^ =- r. Define I{V) by R ■ {£ e R\e • f ^ for all f e V}. If 
J C I{V), J ^ I{V), then we have dimc{x)R/J < dimcV^ = r because of the 
zero-dimensionality of J. Therefore, we have 

J = I{V{J)). 

Under this correspondence of ideals and solutions, a zero-dimensional ideal J 
of R is reducible if and only if there exists a proper subspace W of the solution 
space of V{J) such that < dimc(x)R/I{W) < dimc(x)R/J- In case of one 
variable, the reducibility is equivalent to saying that the generator of the ideal 
can be factored in R. 



Theorem 5.2 The systems of differential equations R ■ Ha{P) is reducible if 
and only z//3 G Z. 

Proof: Any curve is Cohen-Macaulay By applying the theorem of Adolphson 
[0, the holonomic rank of Ha{P) is a„ for all j3. 
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Put M(/3) = An/HA{P)- Consider the left A„-niorphisni 

9i : M{P)-^M{(3+1). (2) 

It has the inverse when /? ^ — 1. Therefore, we have M{—1) ~ M(— 2) ~ 
M(-3) ~ • • • and M(0) ~ Af (1) ~ M{2) ~ • • •. 

When /3 S N, the system admits polynomial solution, then it is reducible. 
It is also easy to see that when (3 S Z<o, the equation is reducible. In fact, 
consider the left Z3-morphism 

di : Af(-l) -^ A/(0). 

It induces a morphism to the solutions by 

/ ^ 5i . /. 

The solution / = 1 of M(0) is sent to zero, so the image of di gives a proper sub- 
space of solutions in the solution space of M(— 1). To find differential equations 
for the subspace, take aU £ such that £di G -H"a(0). Then, {£} C HAi^l). By 
the isomorphism (0), we conclude that when f3 6 Z<oj the system is reducible. 

Let us prove that the system is irreducible when /3 ^ Z by applying the result 
of Beukers, Brownawell, and Heckman Q. For this purpose, we firstly construct 
convergent series solution of Ha(P)- Take w = (1, 1, . . . , 1, 0). Then the degree 
of inyj{lA) is equal to a„. Since I a contains the elements of the form d°'" — d!^', 
the radical of inu,(/A) is {di, . . . , 9„_i). Therefore, the top dimensional standard 
pairs have the form {d^, {n}). 

Let V be the zero of the indicial ideal associated to {d^, {n}): 

Vi = bi,.. . ,W„-1 = bn~l,Vn = . 



Assume f] ^ Z or f3 ^ 0. Taking the lattice basis (a2, — 1,0, . . . ,0), ..., 
(a„, 0, . . . , 0, — 1), we have the following a„ linearly independent convergent se- 
ries solutions 

E 

m6N"- 



f2('^2 -l)---{v2-m2 + l)--- V„{Vn ~ I) ■ ■ ■ {Vn - TTln + 1) 

(wi + 1) (wi + 2) • ■ • {vi + J2 aumk) 

(3) 



X2 J \Xn 

We consider the change of variables: 

Vi = xi,y2 = xl''/x2,. .. ,y„ = x1"/xn- 

The inverse of this change of variables is also rational and the change of variables 
induces that of di. We denote by <& the operation of these change of variables 
of Xi and di. Since the irreducibility is invariant under any birational change of 
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variables, we will prove the irreduciblity of the ideal J = R ■ $(a: '"Ha{(3)x'") 
where R = C{y){dy^,.. ■,dyj. 

Let V be the solution space of J spanned by the series $((p|) • x~") near 
y — 0. If J = I{V) is reducible, then there exists a proper subspace W of 
J such that < dinic(y) R/I{W) < a„. We consider the vector space W = 
{ /(O, . . . , 0, j/„) I / e W}. It is easy to see that dinic(j,„)i?7/(W^') < dimcW^ 
where R' = C{yn){dyJ). Let us prove that dimcW^' = a„ when /9 ^ Z, which 
implies the irreduciblity of J by a contradiction. 

We restrict the series $((||) • x~") to yi = xi = 0, j/2 = x°i 1^2 — 0, . . . , 
y„_i — x-^"^^ /xn — and replace y„ = x^" /xn by x. Without a loss of gener- 
ality, we may assume vi = 0. Then the restricted series has the form 

V" VniVn-l)---iVn-m+l) ^ _^^,„ 

^, ianmy. ^ '> ■ ^^> 

m—O • ' 

It is annihilated by the ordinary differential operator 

{an0z)ianOz - 1) • • • ian9z - a„ + 1) - z{0z + v„). 

By replacing zj a'^ by a;, we obtain the generalized hypergeometric ordinary 
differential equation 

0x{0x ~ l/fln) ■■■[0x~ [Cin - l)/a„) - x(B-j, + U„). (5) 

By B , this ordinary differerential equation of rank a„ is reducible if and only if 
"Vn — k/an ^ Z for all fc = 0, 1, . . . , a„ — 1. If one of them is an integer, /3 becomes 
an integer. Therefore, the ideal I{W') contains the principal ideal generated by 
(||), which is maximal when /3 ^ Z. We conclude that I{W') is generated by (||) 
and hence dime W = a„. [] 
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